Competitive Brownian and Levy walkers 
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Population dynamics of individuals undergoing birth and death and diffusing by short or long 
ranged twodimensional spatial excursions (Gaussian jumps or Levy flights) is studied. Competitive 
interactions are considered in a global case, in which birth and death rates are influenced by all 
individuals in the system, and in a nonlocal but flnite-range case in which interaction affects indi- 
viduals in a neighborhood (we also address the noninteracting case). In the global case one single 
or few-cluster configurations are achieved with the spatial distribution of the bugs tied to the type 
of diffusion. In the Levy case long tails appear for some properties characterizing the shape and 
dynamics of clusters. Under non-local flnite-range interactions periodic patterns appear with peri- 
odicity set by the interaction range. This length acts as a cut-off limiting the influence of the long 
Levy jumps, so that spatial conflgurations under the two types of diffusion become more similar. By 
dividing initially everyone into different families and following their descent it is possible to show 
that mixing of families and their competition is greatly influenced by the spatial dynamics. 

PACS numbers: 05.40.-a, 05.40.Fb, 87.18.Hf, 87.23.Cc 



I. INTRODUCTION 

Birth and death are the most relevant processes in de- 
termining the dynamics of biological populations which 
in the context of statistical physics can be modeled us- 
ing interacting particle models where particle number is 
changing in time. As it is understood by now, birth and 
death processes are also responsible for clustering mech- 
anisms in systems where random- walking individuals un- 
dergo reproduction and death. As a result, aggregation 
of organisms can occur even in simple models where birth 
and death processes are combined with spatial diffusion. 
In fact, in the most simple Brownian bug model, where 
particles reproduce and die with the same probability and 
undergo Brownian motion [l|-y]) clustering of particles 
was observed. In this model the clustering is produced 
simply by the reproductive correlations (the offspring is 
born at the same location of the parent) and by the irre- 
versibility of the death process. Birth and death models 
of moving individuals are the pertinent framework to cap- 
ture properties of biological systems such as planktonic 
populations [^, or patterns in amoebae [4| and bacteria 

Taking into account another central ingredient that is 
present in ecological systems, namely, the competition 
with other individuals in the neighborhood for resources, 
the formation of periodic spatial structures was observed 
in Refs. ^j^. In these nonlocally interacting Brownian 
hug models it was assumed that the reproduction proba- 
bility depends on the number of other organisms in the 
neighborhood. In Ref. [9j nonlocally interacting Levy 
bugs, i.e., reproducing and dying organisms that undergo 
Levy flights, were studied. This type of motion is relevant 
to model cell migration jl(|, biological searching strate- 
gies [ll|, [l2| b acteria dynamics [l^ , or pattern formation 
of mussels |ij| . In Ref. Q it was shown that the forma- 



tion of a periodic pattern is robust with respect to the 
type of spatial motion that the particles perform. The 
periodic arrangement of clusters in these nonlocally in- 
teracting bug models is a consequence of the competitive 
interaction and has a spatial scale determined by the in- 
teraction range |6]. However, a deeper analysis of the dif- 
ferences and similarities between the Brownian and Levy 
cases is still missing. In particular, as shown in [7, 15], 
this analysis can be very conveniently performed by con- 
sidering the limit of the interaction distance reaching the 
system size (global interaction), since a unique cluster 
appears which helps to understand and characterize the 
cluster properties, and the fluctuations of the population 
size. 

In the present paper we report on differences between 
the systems of Brownian and Levy bugs, in the situations 
of global and non-local interactions, as well as in the non- 
interacting case. In addition, results on the dependence 
of population on diffusion, and mixing of families of par- 
ticles are presented for the finite-range interaction case. 
The paper is organized as follows: in Sec. |n]we describe 
the models to be analyzed. In Sec. lIIII thc noninteracting 
bug systems are studied. The infinite competition range 
where each particle is competing with all the others is an- 
alyzed in Sec. IIVI Finally, the nonlocally interacting (i.e. 
with a finite interaction range) models are investigated 
in Sec. Ivl 



II. MODEL AND NUMERICAL ALGORITHM 

We consider a system consisting initially of A^o point- 
like particles, which could be thought as being biological 
organisms or bugs, placed randomly in a two-dimensional 
L X L square domain with periodic boundary condi- 
tions. Except when explicitly stated, we take L = 1, 



so that lengths are measured in units of system size. 
Bugs diffuse, reproduce at rate r^, and die at rate r^; 
i = f , . . . ,iV, and N = N{t) is the number of bugs in 
the system at time t. The numerical algorithm used to 
evolve the system follows the one suggested in Ref. |la |. 
The following sequence of steps is repeated until the final 
simulation time is reached: 

We first compute the random time r after which the 
next demographic event (birth or death) will occur. For 
this we need to determine the total birth and death rates, 
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Btot ^J^^'b, Dtot = Yl '' 



(1) 



i=l 



i=l 



and compute also the total rate 



N 



Rtot = Btot + Dtot ^J2^rl + r^) , 



(2) 



For the random times r we choose an exponential prob- 
ability density with the complementary cumulative dis- 
tribution 



p{t) — exp(— t/t) 



(3) 



so that values of t could be generated from r = — f ln(^o)j 
where ^o is a uniform random number on (0, 1) |17| . The 
characteristic time or time-scale parameter f = (r) is 
determined by the total rate: 



r; 



(4) 



After the random time r, an individual i, chosen 
among all the N{t) bugs, either reproduces or disappears. 
With probability Btot/^tot the event is reproduction and 
with probability -Dtot/^tot it is death. The probability of 
choosing a particular individual i is weighted proportion- 
ally to its contribution to the corresponding total rate. 
In the case of reproduction, the new bug is located at 
the same position {xi,yi) as the parent individual i. Fi- 
nally, all the bugs perform a jump of random length (. in 
a random direction characterized by an angle uniformly 
distributed on (0, 27r) (f and the direction of the jump are 
independent for each particle). The new present time is 
t-fr, bugs are relabeled with indices i ~ 1,2, ..., N{t + T), 
and the process is repeated. 

When bugs undergo normal diffusion (Brownian bugs) , 
a Gaussian jump-length probability density function is 
used. 
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(5) 



with second moment (£^) = £'^; -^ is the space-scale pa- 
rameter. Since we draw the angle specifying the direction 
of the jump from the interval (0, 27r), we restrict i in Eq. 
([5]) to have positive sign. The random jump length £ can 
be computed from £ = £^Gi where S.g is sampled from the 



standard Gaussian distribution with average and stan- 
dard deviation 1, and neglecting the sign. Note that the 
random walk defined in this way is not exactly the same 
as the one in which the walker performs jumps extracted 
from a two-dimensional Gaussian distribution, but it also 
leads to normal diffusion and allows a more direct com- 
parison with the Levy case. The corresponding diffusion 
coefficient can be estimated as 



^={i^)/{2{r))=£y{2f). 



(6) 



As we choose to fix the value of k, and the demographic 
rates, then the space-scale parameter is determined by 



£^ \/2k?= VW^ 



(7) 



In order to simulate the system where the bugs un- 
dergo superdiffusive Levy flights (Levy bugs) one can 
use a symmetric Levy-type probability density function 
for the jump size {£ > 0), behaving asymptotically as 

MM 



Mi)^n£\ 
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^ -> oo {£:$>£) 



(8) 



with the Levy index < /i < 2. For all Levy- type 
probability density functions with /i < 2 the second mo- 
ment diverges, {£'^) = oo, leading to the occurrence of 
extremely long jumps, and typical trajectories are self- 
similar, showing at all scales clusters of shorter jumps 
interspersed with long excursions. For < a < n < 2 
fractional moments {£"') are finite. For the Levy index in 
the range 1 < /i < 2 the value of {£) is finite. The com- 
plementary cumulative distribution corresponding to (|8]) 
behaves as 
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(9) 



As a simple form of complementary cumulative distri- 
bution function which behaves asymptotically as ([9]), we 
use 



P^i£) = {l + b'^''£/£y 



(10) 



with b = [r(l - fi/2)r{p/2)]/T{fi), and £ > 0. As before, 
the direction of the jump is assigned by drawing an ran- 
dom angle on (0,27r). The particular expression for b is 
chosen for consistency with previous work ^] . It gives to 
the tail of the jump distribution the same prefactor as for 
the Levy-stable distribution [20|, but any other positive 
value of b should lead to the same results reported here. 
One can generate a random step-length £ by inverting 
(HUl): 
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with ^0 being a uniform random variable on the unit in- 
terval. Now the diffusion coefficient ^ is infinite, but 
one can define a generalized diffusion coefficient in terms 
of the scales £ and f as 1I8, 191 



K^ = r/(2f ) . 



(12) 



Therefore, in the case of the Levy flights, when fixing the 
value of K^, the space-scale parameter is: 



(2«:,,f)i/^ 
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As we consider the bugs to be point-like, the spatial 
dynamics does not include any interaction between them. 
The interaction is instead taken into account through 
reproduction and death rates, which we assume to be 
affected by competitive interactions. 

If the birth and death rates of a bug are influenced by 
the number of other bugs within a certain radius i?, one 
talks about a nonlocal interaction of finite range. In the 
present paper we assume that the birth and death rates 
of the i-th individual depend linearly on the number of 
neighbors in the interaction range [a|, 



max (O, rbo - 
max (O, rdo 



-liNl 



(14) 
(15) 



Here iV|j. is the number of bugs which are at a distance 
smaller than R from bug i, the parameters rf,o and r^o 
are the zero-neighbor birth and death rates, and the pa- 
rameters a and /3 determine how r^ and r^ depend on 
the neighborhood. For positive values of a and /3, the 
more neighbors an individual has within the radius R, 
the smaller is the probability for reproduction and the 
larger is the probability that the bug does not survive, 
which could arise from competition for resources. The 
function max() enforces the positivity of the rates. Since 
we take R < L/2 (in fact R<^ L), the periodic boundary 
conditions are straightforwardly implemented and bugs 
are never counted twice. 

If the birth and death rates of a bug are instead influ- 
enced by all the other individuals in the system, i.e., 

rl = n, = maii{0, no- a[N{t)-l]} , (16) 
r^ = rd=max{0,rdo + /3[A^(t)-l]} , (17) 

then one talks about global interaction. This is formally 
equivalent to Eqs. p^ and p^ with R sufficiently large 
for the interaction domain to include the whole system, 
but taking care of counting each bug only once, so that 
N}, = N{t) ~ 1. 

In the case the rates of the demographic events are the 
same for all the bugs and assume constant values. 



n = no , 



rd = Tdo , 



P 



(18) 
in Eqs. dH]), ([15l), bugs 



which is equivalent to a 
are noninteracting. 

In the following we discuss the Brownian and Levy bug 
systems when individuals do not influence each other and 
when inter-particle interaction occurs, either global or 
nonlocal. Although we formally maintain the parameter 
/3 in Eqs. (J15l) and (J17p . in our numerical examples we 
restrict to /3 = 0. 



III. SIMPLE BUG MODELS WITH NO 
INTERACTION 

A. Noninteracting Brownian bugs 

The simple Brownian bug model with no interaction, 
i.e., when the birth and death rates of the individuals 
are given by Eq. (|T5| , has been studied and discussed in 
various works [J [2[ . The ensemble average of the total 
population size follows 



{N{t)) = No exp [Mt - to)] 



(19) 



independently of the diffusivity of the bugs; it only de- 
pends on the difference IS. — r\j-- rd- If the birth rate is 
larger than the death rate, A > 0, the total population 
generally explodes exponentially, though there is a finite 
probability for extinction that depends on the initial size 
of the population and decreases with increasing A. If 
the death rate is larger than the birth rate, A < the 
extinction of the population takes place with probabil- 
ity 1. If birth and death are equally probable, A = 0, 
then the average over many realizations is {N(t)) = No 
and the average lifetime is infinite. However, in single 
realizations the fluctuations in the number of individuals 
are huge leading to fast extinction in some runs. In fact, 
there exists a typical lifetime proportional to A^Oi defined 
as the time for which the fluctuations become as large as 
the mean value, beyond which the population is extinct 
with probability close to 1 [1|. 

As a surprising effect, in the systems where the nonin- 
teracting Brownian bugs undergo death and reproduction 
with equal probabilities, spatial clustering of the bugs 
was observed in single realizations [l|-y|- A typical time 
evolution of such a system is illustrated in Fig. [T^. We 
note that in all figures presenting the spatial configura- 
tions of the bugs, we have divided the individuals accord- 
ing to their initial position into nine groups characterized 
by different colors as in Fig. [1] at time t = 0; if reproduc- 
tion takes place, the newborn bug assumes the same color 
as the parent. From Fig. [T^ one can see that many small 
clusters form some time after starting from a uniform 
initial distribution. The occurrence of the clustering is 
related to the fact that in the case of reproduction the 
new bug is located at the same position as the parent. 
Due to the fluctuations and irreversibility of death the 
number of clusters decreases in time, until there will be 
a single cluster consisting of individuals coming from a 
single ancestor. There are constant, spontaneous, short- 
lived break-offs from the main cluster, which are always 
located near it. The center of mass of such a cluster un- 
dergoes a motion similar to that of a single bug 1] and 
its linear width fluctuates with a typical value propor- 
tional to y/No li|. Furthermore, the larger the diffusion 
coefficient k, the wider is the cluster (notice that when 
simulating the system numerically, if the diffusivity be- 
comes so large that the jump lengths become comparable 
to the system size, one needs to take a larger simulation 
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FIG. 1: (Color) Simple bug models with no interaction; spatial configuration of bugs at different times t: (a) Brownian bugs with 
K = 10"'' and (b) Levy bugs with k^ = 10~^ and fi = 1. Reproduction and death occur with equal probability, rj, = Vd = 0.1, 
and the initial number of individuals is A'^o = 1000. Bugs are colored with the color their ancestors bear in the panel at i = 0. 



box). Finally, due to the fluctuations also the last cluster 
disappears. 



IV. GLOBAL INTERACTION 



A. Formation of a cluster 



B. Noninteracting Levy bugs 

In the case of noninteracting Levy bugs, the number 
of individuals still foUovi^s Eq. ([19]) , independently of the 
Levy index /i, and also the cluster formation observed in 
the case of Brownian bugs takes place. Now, however, as 
bugs can perform long jumps, there are also small clus- 
ters continuously appearing and disappearing far from 
the main clusters (Fig. [TJ)). The smaller the value of jjl 
the more anomalous the system, i.e., the larger is the 
probability for long jumps and therefore there are more 
flash-clusters. When the number of clusters has already 
decreased to one, due to the long jumps and fluctuation 
of the number of individuals, new clusters that are placed 
far from the central cluster can appear in the system also 
for some time and often the disappearance of the main 
cluster takes place whereas another new central cluster 
appears somewhere else. As a result the center of mass 
undergoes anomalous diffusion as single bugs do. The 
value of the Levy index /^ influences also the linear size 
of the main clusters: the smaller is /i, the more compact 
are the clusters, although also more particle jumps to 
long distances occur. The influence of the value of k^ 
is similar as in the case of Brownian bugs, i.e., a larger 
value of the anomalous diffusion coefficient results in a 
larger linear size of the clusters. 



Let us now investigate the behavior of the Brownian 
and Levy bug systems in the case of global interaction, 
i.e., birth and death rates of the individuals are given by 
Eqs. (IT5|) . pT|) . The time evolutions of the globally in- 
teracting Brownian and Levy bug systems are illustrated 
by Fig. [2^ and[2jD, respectively. In both systems we start 
from A'o — 500 bugs uniformly distributed in the sim- 
ulation area and choose for the parameters characteriz- 
ing death and birth rates the following values: r^o = 1, 
fdo — 0.1, a = 0.02, /3 = 0. As in the noninteracting 
case, the final state of the dynamics is complete extinc- 
tion, since there is always a nonvanishing probability for 
a fluctuation strong enough to produce that. However, if 
the number of bugs in the system is large this happens 
at very long times [2l|. Then, there is a long-lived qua- 
sistable state for which the average number of individuals 
{N{t)) can be estimated from the condition that death 
and birth are equally probable, r^ — r'^. From there. 
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(20) 



where Aq — rbo — ''do- We have restricted to parameter 
values so that the max functions in Eqs. (116m 7p do not 
operate. Since we have chosen Nq > {N{t)) — 46 in Fig. 
[21 death is more probable at small times and the number 
of bugs decreases rapidly. Approximately at time t = 30 
the number of individuals has reached the value at which 



r = 



r = 20 



r = 50 



/=100 



r = 500 



(a) 



0.5 



(b) 



'^••.••■l>'>r:.-.vv 
i. ••>-C.-'. r'i*"" 


1 • . 


1 
-♦ « 

* 

1 


1 
i 

1 


1 




' ■ . ■ . . " 

■. r 


1 
1 


1 

• 

•• 

1 


1 



0.5 



0.5 



0.5 



0.5 



0.5 



FIG. 2: (Color) Globally interacting bug models; spatial configuration of bugs at different times t: (a) Brownian bugs with 
K — 10"^ and (b) Levy bugs with k^ — 10~^ and /i = 1. The parameters in the reproduction and death rates are: rto = 1, 
Tdo ~ 0.1, a — 0.02, /? = 0. Bugs are colored with the color their ancestors bear in the panel at i = 0. 



death and birth become in average equally probable and 
after this time particle number fluctuates around that 
value; parameters of the birth and the death rates can 
be chosen so that these fluctuations are weak. At this 
time small clusters start to form due to the reproduc- 
tive pair correlations. As in the case of noninteracting 
bugs, fluctuations and irreversibility of death makes the 
number of clusters to decrease in time, although now the 
fluctuations of the particle content of the different clus- 
ters are correlated to keep the total number close to the 
value given by Eq. (|20p and the process is faster. Finally 
a single cluster consisting of bugs coming from the same 
ancestor remains (as stated before, it will also disappear 
at very long times due to finite-size fluctuation effects) 
though there are also now spontaneous short-lived break- 
offs from the central cluster as in the case of noninteract- 
ing bugs. The center of mass of such a cluster is moving 
in space and its linear size is a fluctuating quantity. The 
clustering of the globally interacting bugs was quantita- 
tively discussed in Ref. [7j for the one-dimensional Brow- 
nian bug system. 



B. Fluctuations of the number of bugs 



As indicated by Eq. (j20p . for given values of a and 
/3, the average number of individuals in the system with 
global interaction depends solely on the difference Aq = 
^60 — ''do- It is independent of the concrete values of rbo 
and r^o, as well as of the value of k or k^ and fj,; in fact, it 
does not even depend on whether the system consists of 
Brownian or Levy bugs. Nevertheless, fluctuations of the 
number of bugs do indeed depend on the values of ri,o and 
rdo , even if the difference Aq , and thus the average num- 
ber of bugs, has the same value. To illustrate this, let us 



calculate from the simulations time series the probability 
distribution of the number of individuals in the globally 
interacting Brownian and Levy bug systems. As can be 
seen from Fig. [31 for a given value of Aq, larger values 
of rfco and rdo lead to larger fluctuations. This is a sim- 
ple consequence of the Poisson character of the birth and 
death events for which fluctuations in each of the instan- 
taneous rates are proportional to the mean rates. When 
the distributions are narrow, they are close to Gaussian. 
For larger rates particle number distribution gets broader 
implying that there is an enhanced probability that par- 
ticle number becomes zero at some moment, after which 
bugs become extinct (remember that what is in fact plot 
in Fig. [3] is the numerical particle number distribution in 
the long-lived metastable state before extinction). For 
the present case with Aq = 0.9 and a = 0.02, /3 = 0, 
rate values above the ones shown in Fig. [3] (i.e. r{,o > 2, 
rdo > 1.1) lead to observable extinction after some tenths 
of thousands of steps. An ecological implication of this 
could be the following: one can think of two colonies of 
organisms of the same type, having both the same equi- 
librium size determined for example by the size of the 
living area. Now if in one of the systems the population 
has no enemies and the natural death rate is low, but 
in the other the death rate is higher due to the presence 
of a predator, then the latter system will more proba- 
bly go to extinction sooner due to the presence of larger 
fluctuations. 



The average cluster shape, cluster width, and 
center of mass motion 



Let us keep in the following a = 0.02, /3 = and 
rfco = 1, rdo = 0.1 [the same parameter values as in 
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FIG. 3: (Color online) Probability distribution of the number 
of bugs in globally interacting bug systems. The results are 
numerically obtained from the time series of the particle num- 
ber in the very long-lived state before the fluctuations lead the 
system to the extinction. For all the curves a = 0.02, /3 = 0, 
and the rate difference is Ao = rto — f^o = 0.9, but the rates 
^60, ''do assume different values: (a) r^o = 1, fdo = O.f; (b) 
no = 1.5, rdo = 0.6; (c) rw = 2, rdo = 1-1. The overlapping 
curves correspond to Brownian and Levy bug systems; the 
distributions do not depend on the type of diffusion nor on 
the values of k, k^ or \jl in this globally interacting case. 



Fig. [5] and in Fig. [3] for curve (a)] and study the behavior 
of the cluster formed in the case of a system with global 
interaction defined by Eqs. ([T6 | - p7)) . As mentioned, even 
after the transition from an uniform distribution of bugs 
to a single cluster (and before eventual extinction at large 
times), at some moment the system can consist actually 
of more than one cluster. In such cases we define that all 
the individuals in the system belong to the same cluster, 
even though in the Levy case the distance between the 
bugs (subclusters) can be rather large. In order to avoid 
the boundary effects as much as possible, in Figs. HE the 
linear size of the simulation area was taken as L = 1000 
and to have enough statistics simulations were run until 
t = 5 X 10^. 

Let us start by analyzing the average shape of the clus- 
ter. The average cluster, p{x,y), is obtained setting at 
each time the origin in the center of mass of the cluster 
(distances under the periodic boundary conditions are 
computed under a minimum distance convention) and 
averaging over a long time (after the transition from uni- 
form distribution to one single cluster but before long- 
time extinction) . A one-dimensional cut of it (say across 
X for y = 0, i.e., p(x) = p{x, y — 0)) is shown in Figs. |4] 
and[51 For the case of Brownian bugs, the tail of the aver- 
age cluster is approximately exponential. A pair distribu- 
tion function, which should be related but not identical to 
the average cluster discussed here, was analytically calcu- 
lated in Ref. [la] for a globally interacting Brownian bug 
model of our type but in which total extinction was for- 
bidden. This quantity also displayed a fast decaying tail. 



In the case of Levy bugs the tail of p{x) follows instead 
a power law, p{x) ~ x~^'^^^\ see Fig. HJj, arising from 
the long jumps. Note that, in the present case of circular 
symmetry, the relation p[x,y)dx dy = R{r)(2Tr)~^dr d9 
of p{x, y) with the radial density of the average cluster, 
R{r)^ where r and 9 are the polar coordinates centered 
at the cluster center, implies p{x) = R{r = |a;|)(27r|a;|)~"'^, 
so that the asymptotic behavior of the radial density is 
R{r) ^ r~(^+^\ This is the same asymptotic behavior 
as the individual radial jumps in (|S]) and it is also the 
asymptotic tail of the probability of displacement from 
the original position of nonreprodueing bugs moving by 
Levy flights [la]- We note also that, for k = k^, the 
central part of p{x) is narrower and higher in the Levy 
than in the Brownian bug system, and it is narrower and 
higher the smaller the value of p (see Fig. [5^). This is 
a somehow counterintuitive effect of the Levy motion on 
clusters, already commented in the noninteracting case: 
increasing anomalous diffusion (smaller p) induces larger 
jumps and longer tails, but at small scales it acts as mak- 
ing the cluster more compact. 

The influence of the diffusivity is similar in both sys- 
tems: the larger is the value of n or k^ the more spread 
is the average cluster (see Fig. [5]) . For the Brownian one- 
dimensional case it was shown in Ref. Q that cluster 
width is essentially the distance associated to the Brow- 
nian walk during the lifetime of a bug and their descen- 
dants. Thus, the width increases as k^/^. In the Levy 
case, defining the distance associated to the walk is more 
subtle, since higher moments of displacements diverge. 
But the behavior of lower ones and dimensional analysis 
indicate that typical displacements during a lifetime scale 
as K^ , and then this should determine the width of p{x) 
or R{r) (i.e., the spatial dependence should occur only 



through the combinations [xKf 



-l/Ml 



-1/m 



]). Impos- 



ing additionally that the total number of bugs in the aver- 
age cluster in this global interaction case does not depend 
on particle motion or distribution, and it is thus indepen- 



dent on the value of k^ we have R{r) 
or 



n-^"^F{rK-^^'^), 



p{x) 



2/^ 



G 



,1/^ 



(21) 



with G{u) = F{u)/u. The analogous scaling form for the 
average cluster in the Gaussian diffusion case is 



p{x) 



K^ 1^1/2 j 



(22) 



The insets in Fig. [5] show the validity of these scaling 
forms. 

In addition to the average cluster shape, which gives 
information on the cluster width, it is also interesting to 
study the fluctuations of the cluster width in time. We 
characterize cluster width at each instant of time by the 
standard deviation of the bug positions with respect to 
the center of mass of the cluster at that time. Then, a 
probability density 7r(cr) is constructed from the values 




FIG. 4: (Color online) (a) p{x), the cross-section of the two- 
dimensional particle density of the average cluster in semi- 
log scale; comparison between the Brownian and Levy bug 
systems; k = 10~^, «;^ — 10~^, rto = 1, rao = 0.1, a = 0.02, 
/3 = 0. (b) The tails of p{x) in log-log scale in the case of 
the Levy bug systems for different values of fj.. Solid lines 
correspond to fitting curves oc x~ 



-(2+m) 



of a at different times. Figure IH] shows that in the case 
of globally interacting Brownian bugs the distribution of 
a is short-tailed. In the case of globally interacting Levy 
bugs, in contrast, the distribution for a is characterized 
by tails with a power law decay with exponent — (1 -I- /i). 
This means that in the latter case the cluster width can 
undergo arbitrarily large fluctuations in time. We note 
that the tails in 77(17) decay with the same exponent as the 
radial density R{r) of the average cluster, thus suggest- 
ing that the tails of the average cluster are produced by 
the large fluctuations in the width of the instantaneous 
clusters (which in fact include splitting events). 

The individual motion of bugs drives the behavior of 
the center of mass of the system. Figure [7] depicts the 
probability density p{Acm) of the jump lengths Acm 
perfornred by the center of mass each time the bug mo- 
tion step is executed in the globally interacting Brownian 
and Levy bug systems. For Brownian bugs it is short- 




.50 -0.25 



0.50 



FIG. 5: (Color online) p{x), the cross-section of the two- 
dimensional particle density of the average cluster in semi- 
log scale for different values of diffusivity: (a) Brownian bugs 
and (b) Levy bugs with fi = 1. Other parameters are as in 
Figs. [5] and 21 The insets check the correctness of the scaling 
forms ([22]) (with G{u) = np and u = x/k^^^) and ^^ (with 
G{u) = K^^^p and u = x/zs:^/"). 

tailed. In fact, from the arguments in Ref. [Tf, the cen- 
ter of mass motion of the cluster for globally interacting 
Brownian bugs is characterized by a Brownian process 
with the same diffusion coefficient as the individual bugs. 
In the case of Levy bugs the probability density of the 
jump lengths of the center of mass is described by a dis- 
tribution with a power-law tail with exponent —(1 -I- /it), 
i.e., the center of mass of the cluster formed in the case 
of globally interacting Levy bugs undergoes jumps that 
follow asymptotically the same law as the single bugs, 
Eq. ([5]), and as the radial tails of the average cluster. 
This reflects the fact commented previously that, due to 
the long jumps of the Levy bugs, additional clusters far 
from the main one appear from time to time, strongly 
displacing the center of mass of the system. Due to the 
fluctuations it is even possible that the cluster that used 
to be the main cluster disappears and a new main clus- 
ter forms somewhere else. As a result the center of mass 
motion undergoes the same type of superdifFusion as the 
individual bugs of the system. 

Extending the arguments for the Brownian bugs [7| 
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FIG. 6: (Color online) Probability density 7r(cr) of the stan- 
dard deviation a of the bug positions with respect to the 
center of mass of the cluster in the Brownian and Levy bug 
systems; n = 10""", k^ = 10^^. The distribution is obtained 
averaging over a long time. The curves corresponding to the 
Levy bug systems are well fitted by oc cr~'^^'^' (not shown). 
Other parameters are as in Figs. [J] U andO 
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FIG. 7: (Color online) Probability density p{Acm) of the 
jump lengths of the center of mass in the Brownian and Levy 
bug systems. Same parameter values as in Fig. [6] The tails 
of the curves corresponding to Levy bugs are well fitted by 
oc A^]j^ ** (not shown). 



(which were themselves adapted from the ones in [I'D to 
the Levy case one can heuristically show that the dis- 
tributions of a and Acm sue related. To this aim one 
makes the approximation that the number of bugs in the 
system is constant, say N, instead of being constant on 
average. The center of mass receives a positive contribu- 
tion from the new bugs appearing (at location Xi) due to 
the reproduction between diffusion steps (say at time ti), 
a negative contribution from the bugs disappearing dur- 
ing that time (say from position Xj at time tj), and the 



direct contribution from the Levy jumps £k of all bugs 
present at the diffusion step: 



Acm — 



1 1 1 ^ 



(23) 



ieB 



jeD 



fc=i 



B and D denote the sets of bugs that have been born or 
dead, respectively, between diffusion steps. The two first 
terms can combined in a single one S w ^~^J2p=i^p 
by considering that the two sets have approximately the 
same number of individuals, n. dp = xi — Xj is the dis- 
placement between a pair of these bugs, one just born 
and the other just disappeared, sampled inside the same 
cluster. Then the modulus of each ap should be of the or- 
der of the cluster width a, which fiuctuates in time with 
probability tails ruled by an exponent —(1 + /i). This 
contribution in Eq. ([25)) gives the motion of the center of 
mass due to the birth and death processes. The contri- 
bution from the individual particle jumps is in the last 
term in (|23p , which is a sum of Levy jumps of parameter 
/i so that the tails of the probability density are char- 
acterized by a decay with the same exponent —(1 -I- fi). 
These heuristic arguments imply that the modulus Acai 
will also be distributed with long tails characterized by 
an exponent —(1 + /i), as observed. 



V. NONLOCAL INTERACTION 

A. Formation of a periodic pattern 

In Refs. [g, |^, Q on the nonlocally interacting Brow- 
nian and Levy bugs it was assumed that the birth and 
death rates of the i-th individual are given by Eqs. (|14l) . 
(TTSl) . In the case of Brownian bugs, for small enough 
diffusion coefficient and large enough Aq, the occurrence 
of a periodic pattern consisting of clusters that are ar- 
ranged in a hexagonal lattice was observed (see Fig. [5Ji- 
c) [a, 9 • For large values of the diffusion coefficient such 
periodic pattern is replaced by a more homogeneous dis- 
tribution of bugs (Fig. [5Jl). In the case of Levy bugs, 
since the diffusion coefficient ^ is infinite, one could 
expect that the spatial distribution will not reveal a pe- 
riodic pattern; however, as shown in Rcf. 9], for proper 
parameters periodic cluster arrangements do indeed oc- 
cur (see Fig. [9]). The reason for the divergence of the 
diffusion coefficient in the Levy case is in the statistical 
weight of large jumps. These large jumps have some in- 
fiuence on the characteristics of the pattern formed, but 
the relevant structure is ruled mainly by the interactions 
between individuals. In the Levy bug system however, at 
variance with the Brownian case, even at small values of 
Kp there are many solitary bugs appearing for short time 
periods in the space between the periodically arranged 
clusters due to the large jumps Q, c.f. Figs. IHt- andlH^. 
However, the periodicity of the pattern is of the order 
of R (the interaction range) in both systems, being only 
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FIG. 8: (Color) Interacting Brownian bug model with R — 
0.1, rt,o = 1, Tdo ~ 0.1 and a = 0.02, /? = 0. Spatial configura- 
tion of bugs at time 45000 in systems with different diffusion 
coefficients: (a) k = lO"'^, (b) k = 2 x 10"^, (c) k = 4 x 10"^, 
(d) K — 10"^. The initial configuration of bugs is the same 
as in Figs, [l] and [2] at time t = 0. 



slightly influenced by k or k^ and /x, as demonstrated in 
Refs. [y, ^through a mean-field theory calculation. 

In Ref. Q also the two-dimensional particle density of 
the average cluster, obtained by setting the origin at the 
center of mass of each cluster forming the periodic pat- 
tern and averaging over all the clusters in the simulation 
area and over time, was studied. In both, Brownian and 
Levy bug systems the central part of the average cluster, 
where most of the individuals are concentrated, was well 
fitted by a Gaussian function, but the way the particle 
density decreases when moving away from the center of 
mass of the cluster is rather different. In the Brownian 
case a Gaussian decay provides a good approximation, 
whereas in the Levy case it is close to exponential. The 
comparison with the systems with global interaction, dis- 
cussed in Scc. lIVCl reveals therefore that the interaction 
range R turns the exponential decay into Gaussian and 
the power law decay into exponential. 

For a given value of diffusion coefficient, there exists 
a critical value of Aq below which the system gets ex- 
tinct, independently of a [y]. Above this critical value, 
for every a the increase of Aq results in the increase of 
the average number of bugs, but the pattern formation is 
not much influenced. The latter is, however, true solely 
if Aq increases through the increase of r^o and the death 
rate is low. Namely, as in the case of global interac- 
tion discussed in Sec. lIVBl an increase of the death rate, 
though accompanied by a compensating increase of birth 
rate, leads to larger fluctuations in the particle number. 
In numerical simulations we have observed that the larger 
are the fluctuations in the number of bugs, the more dif- 



FIG. 9: (Color) Interacting Levy bug model with R = 0.1, 
r(,o — 1, Vdo = 0.1 and a = 0.02, P — (same parameters 
as in Fig. [S] for Brownian bugs). The spatial configuration 
of bugs at time 45000 in systems with different generalized 
diffusion coefficients and anomalous exponent: (a) fi:,i = 10^*, 
/i = 1; (b) K^ = 10-^ M = 1; (c) K^ = 10-^ fi = 1.5; (d) 



5 X 10" 



1.5. The initial configuration of bugs is 



the same as in Figs. [T] and [2] at time t — 0. 



flcult is the formation of the periodic pattern, and flnally 
the individuals do not arrange in the periodic pattern but 
in random clusters (see also Ref. [la|)- This effect may 
in fact make difficult to observe the periodic clustering 
phenomenon in real competitive biological systems. 

In the following we keep for the parameters in the birth 
and death rate the same values as in the case of global 
interaction, i.e., rf,o = 1, rdo = 0.1, a = 0.02, P — 0. 
For these parameter values the number of bugs fluctu- 
ates only weakly around the mean value. Differently from 
the case of global interaction, now the average number of 
bugs in the system is influenced not only by the birth and 
death rates, but also by the diffusion, i.e., in the case of 
Brownian bugs by k and in the case of Levy bugs by k^ 
and n, see Fig. (TUl The smaller is k, k^ or /i, the larger 
the particle number. At the same time Figs. [8] and [9] in- 
dicate that by decreasing k or k^ the linear width of the 
clusters becomes smaller, the particle density in the clus- 
ters higher, and the density between the clusters lower 
(c.f. Sec. II V Cl and see also Ref. 9]). Somehow counter- 
intuitively, the effect of decreasing ^ seems to have the 
same effects, as commented above for the noninteracting 
and global cases. Furthermore, the value of k or k^ and 
/i seems to weakly influence the number of clusters in the 
system: In Figs. and[8]and[9]smaller values lead to larger 
number of clusters. This observation is not explained by 
the linear instability analysis of [g, Q . 
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FIG. 10: (Color online) a): Average number of bugs versus 
diffusion coefficient in the system with Gaussian jumps, b): 
Average number of bugs versus anomalous exponent ^ in the 
system with Levy jumps for various values of the anomalous 
diffusion coefficient. Other parameters as in Figs. [8]and[9l 



B. Mixing of different families 

It is interesting to study the evolution of the system 
also regarding the disappearance or survival of the dif- 
ferent groups, by dividing initially the bugs into differ- 
ent families and following their descent. In the case of 
nonlocally interacting Brovifnian bugs, a very low diffu- 
sion coefficient leads to the situation in which after clus- 
ter formation the inter-cluster travel is very rare because 
the individuals are not capable to make the jumps from 
one cluster to another one, and it is also very unlikely 
to arrive to the next cluster doing a multistep random 
walk because death is very probable between the clus- 
ters. Therefore, in the case of very low diffusion different 
families would remain inside their initial clusters. If one 
assumes that initially each individual represents a dif- 
ferent family, then only inter-cluster competition occurs 
and the final number of families is equal to the number 
of clusters. If instead initially individuals are assigned 



to families according to large areas of initial positions 
(larger than typical cluster size as done in Figs. [T] and 
[2] at time i = 0) , there is no family competition inter- 
nal to the clusters, most families survive and the clusters 
coming from different families occupy approximately the 
territory of the ancestors even after a long time, as can be 
seen from Fig. [SK- In that case, the travel of a cluster to 
a new territory away from the other clusters of the same 
family can take place due to the diffusion of the cluster 
as a whole during the clusters arrangement into the pe- 
riodic pattern. For larger values of k the inter-cluster 
travel is possible which leads to the conquering of new 
territories, i.e., bugs can be found in a region where their 
ancestors were not from, Fig. ^p. The effect is larger for 
larger k and leads to the disappearance of some families, 
as can be seen from Fig. [St. Finally, for increased dif- 
fusion, intra-cluster competition will force all surviving 
bugs to be from a single family (in fact, from a single an- 
cestor); which family (ancestor) wins is a random event. 
The process is faster for larger diffusion. Increasing the 
diffusivity further even the periodic pattern disappears, 
Fig.[5H. 

Figure [S] illustrates the family mixing for nonlocally in- 
teracting Levy bugs. In this case the inter-cluster travel- 
ing is supported by the long jumps. Differently from the 
case of Brownian bugs, now the individuals can reach 
not only the next neighboring cluster but also clusters 
far away. Consequently, a cluster originally consisting of 
bugs coming from one ancestor can after some time turn 
into a cluster consisting of bugs coming from different 
families placed initially far away. Thus, intra-cluster bug 
competition becomes soon competition between families, 
and even if the diffusivity of the bugs is very low, at the 
end the Levy bug system would consist of individuals 
coming from one or just a few ancestors. As in the case 
of Brownian bugs, the process of the disappearance of 
families is faster the greater is the generalized diffusion 
coefficient. 

Besides the diffusion of a cluster as a whole during the 
formation of the periodic pattern and the conquering of 
new territories through the migration to and survival in 
another cluster, the mixing of clusters from different fam- 
ilies can take place also due to the appearance of a new 
cluster if in the periodic pattern there is a dislocation. 
In the case of Brownian bugs the new cluster is formed 
through the splitting of an old cluster. In the case of 
Levy bugs, instead, the new cluster can appear also far 
from the original territory. 



VI. CONCLUSIONS AND OUTLOOK 

We have presented some detailed properties of interact- 
ing particle systems in which the individuals are Brown- 
ian or Levy random walkers which interact in a compet- 
itive manner. We have seen strong differences between 
the globally and the finite-range nonlocally interacting 
systems. In the systems with global interaction the spa- 
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tial distribution of the bugs becomes tied to the type of 
diffusion, Brownian or Levy. Typical configurations con- 
sist of a single or a few clusters for both types of motion. 
For the Levy bug systems long tails appear in the mean 
cluster shape and in probability distributions of cluster 
width and of jumps of the center of mass. For Brownian 
bug systems these quantities appear to be much shorter 
ranged. This is qualitatively also the situation in the non- 
interacting case, although then the effects of the particle- 
number fluctuations are much stronger. Under non-local 
finite-range interactions the situation is rather different. 
First, single cluster configurations are generally replaced 
by periodic patterns with periodicity set by the interac- 
tion range R. Motion of individual clusters is severely 
restricted by the presence of the neighboring clusters. In 
addition, the natural spatial cut-off introduced by the in- 
teraction range R seems to limit the influence of the long 
Levy jumps, so that measures of spatial cluster shape 
do not generally exhibit power laws, making spatial con- 
figurations under both types of diffusion more similar. 
Mixing of families and their competition is nevertheless 
greatly influenced by the type of motion. This suggests 
that it would be interesting to consider the influence of 
different types of diffusion into competitive genetic mix- 
ing processes |22| . 



Obtaining analytic understanding in this type of inter- 
acting systems is difficult, but at least the nature of the 
instability leading to pattern formation and its relevant 
spatial scale have been clarified by using partial integro- 
differential equation descriptions of the mean field type 
[g, Q , which are useful in broader contexts |23l 4 2a l . How- 
ever, from previous work in the Gaussian case |7|, |8|, |26l| , it 
is known that quantities such as cluster width and struc- 
ture or transition thresholds strongly depend on particle- 
number fluctuations. Thus, obtaining additional results 
from differential equation approaches would need the in- 
clusion of effective multiplicative noise terms [a] or focus 
on statistical quantities such as pair correlation functions 
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